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We derive the Boltzmann equation in the synchronous gauge for massive neutrinos with a deformed 
dispersion relation. Combining the 7-year WMAP data with lower-redshift measurements of the 
expansion rate, we give constraints on the deformation parameter and find that the deformation 
parameter is strong degenerate with the physical dark matter density rather than the neutrino mass. 
Our results show that there is no evidence for Lorentz invariant violation in the neutrino sector. 
The ongoing Planck experiment could provide improved constraints on the deformation parameter. 
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I. INTRODUCTION 


Neutrino oscillations imply that at least two of the 
three neutrino types have non-zero mass (see Ref. [1] 
and reference therein). Unfortunately, neutrino oscilla- 
tion experiments only provide two mass squared differ- 
ences, but not the overall mass scale. Cosmology provides 
a promising way to tackle this problem by the gravita- 
tional effect of massive neutrinos on the expansion his- 
tory near the epoch of matter-radiation equality and on 
the growth of large-scale structures at late time. Preci- 
sion measurements of the cosmic microwave background 
(CMB) anisotropies and the large-scale structure distri- 
bution of matter have allowed us to determine or con- 
strain the absolute neutrino mass scale [2-5] (see also 
Ref. [6] for a recent review). 


Indeed, neutrino oscillations can be explained by small 
Lorentz invariance violation even without introducing 
neutrino mass, as shown in [7, 8, 10-12]. Observed neu- 
trino oscillations may originate from a combination of 
effects involving neutrino masses and Lorentz invariance 
violation. The possibilities of Lorentz invariance viola- 
tion were explored in quantum gravity [13], loop quan- 
tum gravity [14], non-commutative field theory [15], and 
doubly special relativity theory [16]. 

It is well-known that Lorentz symmetry is a funda- 
mental feature of modern descriptions of the nature, in- 
cluding both the Einstein’s general relativity and stan- 
dard model of particle physics. One can expect that 
the breaking of Lorentz symmetry may leave some im- 
prints in astrophysical observations such as the CMB 
anisotropies and large-scale structure of our universe. In 
this paper, we consider cosmological tests of Lorentz in- 
variance violence in the neutrino sector. Here we focus 
on the Coleman-Glashow model, in which the energy- 
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momentum relation is modified by a Lorentz-violating 
interaction in the framework of conventional quantum 
field theory [8]. We will build up a Lagrangian and de- 
rive the Boltzman equation in the synchronous gauge for 
the massive neutrinos with deformed dispersion relation, 
and then place constraints on Lorentz invariance viola- 
tion by combining the 7-year WMAP data [3] with the 
latest distance measurements from the Baryon Acoustic 
Oscillation (BAO) in the distribution of galaxies [17] and 
the Hubble constant (Ho) measurement [18]. 

This paper is organized as follows. In Sec. II we write 
down the Lagrangian for neutrinos with deformed disper- 
sion relation. In Sec. III we derive the Boltzmann equa- 
tion for neutrinos in the synchronous gauge. In Sec. IV 
we place constraints on the deformation parameters us- 
ing the CMB data in combination with measurements of 
the Hubble constant Ho and the BAO feature. Section V 
is devoted to conclusions. 


II. DEFORMED DISPERSION RELATION 


At a phenomenological level the deformed dispersion 
relation for massive neutrinos can be generally parame- 
terized by 


m 


E? =m? +p +Y ang: (1) 
n=1 


where FE is the neutrino energy, m is the neutrino mass, 
p = (p'p;)\/? the magnitude of the 3-momentum, a&n’s are 
dimensionless coefficients and M denotes the energy scale 
corresponding to Lorentz symmetry violation (which is 
typically taken to be the Planck mass). Such a deformed 
dispersion relation implies that there are departures from 
Lorentz invariance in the neutrino sector if a, Æ 0. The 
n = 1 term would produce huge effects at low energy and 
has been strongly constrained. The p” term with n > 3 
is suppressed by 1/M”~?. In the present work we will 
therefore consider the case of n = 2, 


E? =m? +p’ +€p’, (2) 
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where € is the deformation parameter characterizing the 
size of Lorentz symmetry violation. Such a deformed 
dispersion relation was constructed by Coleman and 
Glashow in the framework of conventional quantum field 
theory [8]. 

Here we point out that the dispersion relation given in 
(1) and (2) is not very general. It neglects oscillations, 
possible species dependence, anisotropies associated with 
violations of rotation symmetry, and CPT violation. As 
shown recently by Kostelecky and Mewes, all of these are 
possible [9]. For example, it is known that taking odd 
values of n in Eq. (1) corresponds to CPT violation and 
gives rise to a sign difference for neutrinos and antineutri- 
nos [9]. We have to emphasize that the model considered 
in this paper is one of many possible Lorentz-violating 
theories. 

The scalar perturbations of the Friedmann-Lemaitre- 
Robertson-Walker metric in the synchronous gauge can 
be written as 


ds” = a*(r) [—dr? + (dij + hij)dx' dx] , (3) 


where a(r) is the scale factor, 7 is the conformal time, 
and the scalar mode of h;; is represented by two functions 
h and 7 which are defined by 


@Bke** Goa T) + (kik; 
1 > 
—55i) 6n(k,7)|, (4) 


in Fourier space, where k = kk. The action for a neutrino 
with the dispersion relation (2) can be written in the first 
order formalism by 


S= fo [2" P, — À (m? + g” P,P, + Eg” P;P;)] , (5) 
where À is a Lagrange multiplier which enforces the mass- 
shell condition (2) and 

Po m —aE, (6) 
1 
P, = aig + shiz)’. (7) 
The above action implies that the Hamiltonian van- 
ishes. Performing the Legendre transformation to the 


Lagrangian formalism and eliminating À from the action 
(5), we obtain the Lagrangian for neutrino as follows 


_ i u? hijuřu 


where u’ = da'/dr is the coordinate velocity and u? = 
d;ju'u?. Following [19] from the Lagrangian (8), the 
equations of motion in the synchronous gauge become 


dx’ O (1+ £)6%q; 


i, 7 [1+O(A)] , (9) 
d i 4 i 
= = Sgin n, (10) 


where q; = qn; is the comoving 3-momentum written 
in terms of its magnitude and direction with n'n; = 


iznni = 1, and € = \/m?a? + (1+ €)q? is the comoving 


energy of neutrinos. 


Ill. BOLTZMANN EQUATION 


The number density n,, energy density p, and pressure 
P, for massive neutrinos with (2) are given by 


3 

m = of Gala) (11) 
3 

m= af Erho. (12) 
3 2 

p = oof no a 


Here the zeroth-order phase space distribution is well ap- 
proximated by the relativistic Fermi-Dirac distribution 
fola) = gs|[1 + exp( V1 + €¢/To)|~*, where Tp is the neu- 
trino temperature today and g, = 2 the number of spin 
degrees of freedom. If m > To, the total mass of neutri- 
nos J m = 94(1+ €)3/20,h? eV. 

The perturbed energy density, pressure, energy flux, 
and shear stress in the Fourier space k are respectively 
given by 


ia = a | Ehon (14) 
ip = a | E HEE now 09) 
oot PO = A f EVT (06) 
(vt Prov a AH OF OH pola), (17) 


3a4 J (27)3 € 


where the perturbations Y; satisfy the following Boltz- 
mann equation 


. q 
P+ (lp D AFI (2 + 1)Wi41 — lY-1] 
fe. DB, «le \ ding 
691 — =n — dqr= = 
+ (Sarit one surgi) ding 0 , (18) 


in the synchronous gauge. Such a Boltzmann hierar- 
chy is effectively truncated by adopting the following 
scheme [20] 


(2lmax + 1)€ 


Ww yy a —wW Zie 19 
Imax t1 © T y gyger Yiee T Pimat (19) 


The initial conditions for the perturbation Y; and 7 on 
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super-horizon scales (kr < 1) are given by 
1. dln fo 


Wo = ~qov ding’ (20) 
€ dln fo 
UV = -— I, — , 21 
: 3/1 +éqk ” ding 4 
7 1 1.. dl fo 
M2 = “oe a i 
5+9V/I FER, — 5R, 7 
= 20 — IV Tsw ee, (23 
" sastavi FER) OO 9) 
where 
ip = -2O(kr), (24) 
2(2+ Ry — VI FER, ) 2 
p = LOOO Olr’, 25 
. 3(1b+4/TFER,) ©) ve) 
(23+ 4R )V1 +£ 


4,3 
w= sedi N S 
Here, C is a dimensionless constant determined by the 
amplitude of the fluctuations from inflation and Ry, = 
pv/(pv + py) during radiation domination. 

In order to compute the theoretical CMB power spec- 
trum, we modified the Boltzmann CAMB code in [27] 
to appropriately incorporate the Lorentz-violating term 
in the neutrino sector. Actually the Lorentz-violating 
term affects not only the evolution of the cosmological 
background but also the behavior of the neutrino pertur- 
bations. From Eqs. (11)-(13), we see that increasing £ 
decreases the number density, energy density and pres- 
sure of neutrinos, and thereby increases the redshift when 
the matter density equals the radiation density and re- 
duces the expansion rate prior to and during the epoch of 
photon-baryon decoupling. This leads to reduced heights 
of the first and second peaks of the CMB while a nearly 
constant increase in the acoustic oscillation amplitudes 
at l > 600. On the other hand, the coefficient of the sec- 
ond term in the Boltzmann equation (18) plays an active 
role, which alters the shape of the CMB power spectrum 
caused by changing neutrino propagation. Decreasing € 
increases fluctuation power both at l < 10 and at l > 100. 
Moreover, the CMB is more sensitive to negative values 
of € than positive ones. These two effects can be distin- 
guished from a change in the total mass of neutrinos or in 
the effective number of extra relativistic species [21-26], 
as shown in Figure 1. 


IV. COSMOLOGICAL CONSTRAINTS 


In our analysis we use a modified version of the pub- 
licly available CosmoMC package to explore the param- 
eter space by means of Monte Carlo Markov chains tech- 
nique [28]. We consider a flat ACDM models plus three 
Lorentz-violating neutrino species, described by a set of 
cosmological parameters 


{Qyh?, Qh’, Os, T, Ms, As, um, E}, 
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FIG. 1: Theoretical angular power spectrum of the CMB for 
€ = —0.1, 0, 0.1. Here X m = 3 x 0.3 eV is fixed. 


where h is the dimensionless Hubble parameter such that 
Ho = 100h kms~! Mpc™, 05h? and Qh? are the physi- 
cal baryon and dark matter densities relative to the crit- 
ical density, ©, is the ratio of the sound horizon to the 
angular diameter distance at the photon decoupling, T is 
the reionization optical depth, n, and A, are the spectral 
index and amplitude of the primordial curvature pertur- 
bations at the pivot scale kọ = 0.002 Mpc™t, X m is 
the total mass of neutrinos assuming that the three neu- 
trino masses are approximately degenerate, and € is the 
deformation parameter. 

We use the seven-year WMAP (WMAP?) data with 
the likelihood code supplied by the WMAP team. We 
consider the Sunyaev-Zel’dovich (SZ) effect, in which 
CMB photons scatter off hot electrons in galaxy clusters. 
Given a SZ template, the effect is described by a SZ tem- 
plate amplitude Agz as in the WMAP paper [3]. We also 
impose Gaussian priors on the Hubble constant, Hp = 
74.2 + 3.6 kms~! Mpc~!, measured from the magnitude- 
redshift relation of low-z type Ia supernovae [18], and on 
the distance ratios, r;/Dy(z = 0.2) = 0.1905 + 0.0061 
and rs/Dy(z = 0.35) = 0.1097 + 0.0036, measured from 
BAO in the distribution of galaxies [17]. Here rs is the 
comoving sound horizon size at the baryon drag epoch 
and Dy is the effective distance measure for angular di- 
ameter distance. 

In Table I we summarize our results. With the data 
of WMAP7+Ho+BAO, the deformation parameter is es- 
timated to be € = —0.077 + 0.089, which implies a null 
detection of Lorentz invariant violence within the error 
limits. Such large uncertainties in € mainly come from 
the strong correlation between the deformation param- 
eter and the physical dark matter density as shown in 
Figure 2. In this case the uncertainties in Qh? are about 
three times larger than those derived in the standard 
ACDM model [3]. Since the deformation parameter is 
nearly uncorrelated with the total mass of neutrinos in 


201709.00166v1 


chinaXiv 


parameter WMAP+Ho+BAO WMAP+Ho+BAO  Planck+Ho+BAO 
1009h? 2.312 + 0.078 2.314 + 0.081 2.260 + 0.023 

Qh? 0.1200 + 0.0103 0.1206 + 0.0104 0.1133 + 0.0033 

Qa 0.715 + 0.020 0.723 + 0.017 0.725 + 0.005 

Ho 70.9 + 2.4 km/s/Mpc 72.1 + 2.2 km/s/Mpc 70.3 + 1.0 km/s/Mpc 
T 0.089 + 0.014 0.088 + 0.014 0.088 + 0.004 

Som < 0.51 eV 0 (fixed) 0 (fixed) 

£ —0.077 + 0.089 —0.073 + 0.089 —0.005 + 0.037 


TABLE I: Mean values and marginalized 68% confidence level for the deformation parameter and other cosmological parameters. 


For the total mass of neutrinos, the 95% upper limit is given. 


Figure 2, our constraints on the deformation parameter 
are not significantly changed if three neutrinos are mass- 
less, as we can see from Table I. Compared to parti- 
cle physics experiments, cosmological observations yield 
much weaker constraints on the Lorentz-violation param- 
eter in the neutrino sector. As listed in Table XIV of 
Ref. [29], previous constraints range from parts in 10° 
to parts in 10! from time-of-flight measurements and 
various threshold analyses. 

Moreover, we present the constraints on the deforma- 
tion parameters from the ongoing Planck experiment [30] 
in Table I. Following the MCMC method described in 
Ref. [31], we generate synthetic data for the Planck ex- 
periment and then perform a systematic analysis on the 
simulated data. As we can see in Table I, the Planck data 
plus measurements of the Hubble constant and the an- 
gular diameter distance will reduce the uncertainties in € 
by a factor of 2.4. Therefore, Planck CMB measurements 
allow us to detect the signature of Lorentz invariant vio- 
lence at 2ø confidence level if || > 0.074. 


V. CONCLUSIONS 


We studied the cosmological consequences of Lorentz- 
violating neutrinos. We obtained the generalized La- 
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FIG. 2: Two-dimensional joint marginalized constraints (68% and 95% confidence level) on the deformation parameter €, 
physical dark meter density Qh? (left) and the total mass of neutrinos $ m (right), derived from the WMAP7+H +BAO 
data. The dashed line corresponds to Lorentz invariance. 
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